Decoherence and asymptotic entanglement 
in open quantum dynamics 



Aurelian Isar 
Department of Theoretical Physics 
National Institute of Physics and Nuclear Engineering 
Bucharest- Magurele, Romania 
e-mail: isar@theory.nipne.ro 

Abstract 

In the framework of the theory of open systems based on completely positive 
quantum dynamical semigroups, we determine the degree of quantum decoher- 
ence of a harmonic oscillator interacting with a thermal bath. It is found that 
the system manifests a quantum decoherence which is more and more signif- 
icant in time. We also calculate the decoherence time and show that it has 
the same scale as the time after which thermal fluctuations become comparable 
with quantum fluctuations. Then we solve the master equation for two inde- 
pendent harmonic oscillators interacting with an environment in the asymptotic 
long-time regime. We give a description of the continuous-variable asymptotic 
entanglement in terms of the covariance matrix of the quantum states of the con- 
sidered system for an arbitrary Gaussian input state. Using the Peres-Simon 
necessary and sufficient condition for separability of two- mode Gaussian states, 
we show that the two non-interacting systems immersed in a common environ- 
ment become asymptotically entangled for certain environments, so that in the 
long-time regime they manifest non-local quantum correlations. 

Key words: Open systems, quantum decoherence, quantum entanglement, insep- 
arable states. 

1 Introduction 



By quantum decoherence (QD) [TJ [2] we understand the irreversible, uncontrollable 
and persistent formation of quantum correlations (entanglement) of a system with its 
environment [3], expressed by the damping of the coherences present in the quantum 
state of the system, when the off-diagonal elements of the density matrix decay below 
a certain level, so that this density matrix becomes approximately diagonal. 

Since QD strongly depends on the interaction between the system and its external 
environment [HE], its role became relevant in many interesting physical problems. 
In many cases one is interested in understanding QD because one wants to prevent 
decoherence from damaging quantum states and to protect the information stored in 



quantum states from the degrading effect of the interaction with the environment. 
QD is also responsible for washing out the quantum interference effects which are 
desirable to be seen as signals in experiments and it has a negative influence on many 
areas relying upon quantum coherence effects, in particular QD is a major problem in 
quantum optics and physics of quantum information and computation [3]. 

Quantum entanglement represents the physical resource in quantum information 
science which is indispensable for the description and performance of such tasks like 
teleportation, superdense coding, quantum cryptography and quantum computation 
pE]. Therefore the generation, detection and manipulation of the entanglement contin- 
ues to be presently a problem of intense investigation. 

When two systems are immersed in an environment, then, besides and at the same 
time with the QD phenomenon, the external environment can also generate a quantum 
entanglement of the two systems and therefore an additional mechanism to correlate 
them [3 El [7] . In certain circumstances, the environment enhances entanglement and 
in others it suppresses the entanglement and the state describing the two systems 
becomes separable. The structure and properties of the environment may be such 
that not only the two systems become entangled, but also such that a certain amount 
of entanglement survives in the asymptotic long-time regime. The reason is that even 
if not directly coupled, the two systems immersed in the same environment can interact 
through the environment itself and it depends on how strong this indirect interaction 
is with respect to the QD, whether entanglement can be generated at the beginning 
of the evolution and, in the case of an affirmative answer, if it can be maintained for 
a definite time or it survives indefinitely in time [5]. 

In this work we study QD of a harmonic oscillator interacting with an environ- 
ment, in particular with a thermal bath, in the framework of the theory of open 
quantum systems based on completely positive dynamical semigroups. We determine 
the degree of QD [8] for different regimes of the temperature of environment. It is 
found that the system manifests a QD which in general increases with time and tem- 
perature. We also calculate the decoherence time and show that it has the same 
scale as the time after which thermal fluctuations become comparable with quantum 
fluctuations. Then we investigate, in the same framework, the existence of the con- 
tinuous variable asymptotic entanglement for a subsystem composed of two identical 
harmonic oscillators interacting with an environment. We are interested in discussing 
the correlation effect of the environment, therefore we assume that the two systems 
are independent, i.e. they do not interact directly. The initial state of the subsystem 
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is taken of Gaussian form and the evolution under the quantum dynamical semigroup 
assures the preservation in time of the Gaussian form of the state. We only investi- 
gate here the asymptotic behaviour of the subsystem states. The time evolution of 
the entanglement, in particular the possibility of the so-called "entanglement sudden 
death", that is suppression of the entanglement at a certain finite moment of time, 
will be discussed in a future work. 

The organizing of the paper is as follows. In Sec. 2 we review the Markovian 
master equation for the damped harmonic oscillator and solve it in coordinate repre- 
sentation. Then in Sec. 3 we investigate QD and calculate the decoherence time of 
the system. In Sec. 4 we write the equations of motion in the Heisenberg picture for 
two independent harmonic oscillators interacting with a general environment. With 
these equations we derive in Sec. 5 the asymptotic values of the variances and covari- 
ances of the coordinates and momenta which enter the asymptotic covariance matrix. 
Then, by using the Peres-Simon necessary and sufficient condition for separability of 
two- mode Gaussian states [9j [10], we investigate the behaviour of the environment 
induced entanglement in the limit of long times. We show that for certain classes of 
environments the initial state evolves asymptotically to an equilibrium state which is 
entangled, while for other values of the parameters describing the environment, the 
entanglement is suppressed and the asymptotic state is separable. The existence of 
the quantum correlations between the two systems in the asymptotic long-time regime 
is the result of the competition between entanglement and QD. A summary is given 
in Sec. 6. 



2 Markovian master equation for a harmonic oscil- 



In the axiomatic formalism based on completely positive quantum dynamical semi- 
groups, the irreversible time evolution of an open system is described by the following 
general quantum Markovian master equation for the density operator p(t) [11] : 



The harmonic oscillator Hamiltonian H is chosen of the general quadratic form 



lator 



dp{t) 
dt 



(1) 




(2) 
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and the operators Vj, Vj , which model the environment, are taken as linear polynomials 
in coordinate x and momentum p x . Then the master equation (OQ) takes the following 
form [T2] : 

dp it % 

-r t = ~j\ H ^ P\ - ^"( A + V)l x ' PVx + Pxp] + ^(A - n)\p x , px + xp) 

Dr ' )x [x, [x,p]} - ^[ Px , [ Px ,p]} + [p x ,p]] + [p x , [x,p]}). (3) 



The diffusion coefficients D xx , D PxPx , D xPx and the dissipation constant A satisfy the 
fundamental constraints: D xx > 0, D PxPx > and 

D xx D PxPx — D 2 xpx > — ^— • (4) 

In the particular case when the asymptotic state is a Gibbs state pc(oo) = e~~kr /Tre~ St , 

these coefficients become 

X — ah Til) X + u too 

D xx = — coth — — , D PxPx = ——hmuj coth — — , D xpx = 0, 5 

2 muj 2k L 2 2k 1 

where T is the temperature of the thermal bath. In this case, the fundamental con- 
straints are satisfied only if A > p and 

(A 2 -/i 2 )coth 2 |^>A 2 (6) 

and the asymptotic values a xx (oo), cr PxPx (oo), a xPx (oo) of the dispersion (variance), 
respectively correlation (covariance), of the coordinate and momentum, reduce to [12] 

/ \ ft . Tiuj Timuj %u , . 

axx ^°°> = 2jfeT' "i"*^ 00 ) = ~~ 2~~ coth 2fcT' ^C 00 ) = °- ( 7 ) 

We consider a harmonic oscillator with an initial Gaussian wave function (<r x (0) 
and <J Px {0) are the initial averaged position and momentum of the wave packet) 

representing a correlated coherent state [13] (squeezed coherent state) with the vari- 
ances and covariance of coordinate and momentum 

hS hmu hr 
°"xx(0) — t^ZTT' a PxpM = —^ —i a x P M = ~ n ~ ( 9 ) 



2mu PxPxy J 25(1 -r 2 )' ^ > 2^\ 



r 



Here S is the squeezing parameter which measures the spread in the initial Gaussian 
packet and r, with |r | < 1 is the correlation coefficient. The initial values (El) correspond 
to a minimum uncertainty state, since they fulfil the generalized uncertainty relation 
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Mo)<w(o) -<J°) = ip ( 10 ) 



For 5 = 1 and r = the correlated coherent state becomes a Glauber coherent state. 
From Eq. (J3j) we derive the evolution equation in coordinate representation: 

dp %h . d 2 d 2 . imuj 2 . 2 , 2 \ 

-^(x - x'fp + D xx (^- + A) 2p _ 2iD xpx %{x - x')(^ + (11) 

The first two terms on the right-hand side of this equation generate the usual Liou- 
villian unitary evolution. The third and forth terms are the dissipative terms and 
have a damping effect (exchange of energy with environment). The last three are 
noise (diffusive) terms and produce fluctuation effects in the evolution of the system. 
D PxPx promotes diffusion in momentum and generates decoherence in coordinate x - 
it reduces the off-diagonal terms, responsible for correlations between spatially sepa- 
rated pieces of the wave packet. Similarly D xx promotes diffusion in coordinate and 
generates decoherence in momentum p x . The D xPx term is the so-called "anomalous 
diffusion" term and it does not generate decoherence. 

The density matrix solution of Eq. ( flTl) has the general Gaussian form 

< x\p(t)\x' >= ( — -r-r)5 exp[ ^— -( X + X - a x (t)) 2 

2ixo xx (t) 2cr xx (t) 2 

~^%-^ X ~ X '? + iH^^T" - a x {t)){x - x>) + U x (t)(x - x% (12) 
2h a xx (t) na xx (t) 2 h 

where 

a(t) = a xx (t)a PxPx (t) - a 2 (t) (13) 



is the determinant of the covariance matrix 

<Txx{t) °xpM 

<?x Px {t) cr PxPx {t) 



(14) 



and represents also the Schrodinger generalized uncertainty function. In the case of 
a thermal bath we obtain the following stationary state solution for t — > oo (e = 
Tiu/2kT): 

TflUJ i 7TIUJ (cC ~\~ x'l^ 

< x]pioo)lx ' >= exp{ - Is'^thT + (I " *'> 2 » th * < 15 > 
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3 Quantum decoherence 



An isolated system has an unitary evolution and the coherence of the state is not lost - 
pure states evolve in time only to pure states. The QD phenomenon, that is the loss of 
coherence or the destruction of off-diagonal elements representing coherences between 
quantum states in the density matrix, can be achieved by introducing an interaction 
between the system and environment: an initial pure state with a density matrix which 
contains nonzero off-diagonal terms can non-unitarily evolve into a final mixed state 
with a diagonal density matrix. 

In the literature several measures of the degree of decoherence have been intro- 
duced, expressed in terms of the von Neumann entropy [H] , linear entropy [T5l [16] or 
the coefficients of the matrix elements of the statistical operator with respect to the 
energy basis [TT] . 

Using new variables £ = (x + x')/2 and A = x — x', the density matrix ([121) 
becomes 

p(S, A, t) = ,/^exp[-a£ 2 - 7 A 2 + t(3ZA + 2aa x {t)Y, + - (3a x {t))A - mx 2 (t)],(16) 

V 7r a 

with the abbreviations 

a = 1 7= a ^ (3 = , (17) 

2o xx {ty 2h 2 a xx (tY Tia xx {t) 

In the present work we use the representation-independent measure of the degree 
of QD [8], which is given by the ratio of the dispersion l/y/2j of the off-diagonal 



element p(0, A, t) to the dispersion y2/a of the diagonal element /?(£, 0, t) 

....... h 

S QD (t) 



It can easily be shown that 6qd is related to the linear entropy [18J. 

The finite temperature Schrodinger generalized uncertainty function has the ex- 
pression [19] (with the notation fl 2 = u 2 — /i 2 , u> > ji) 

a(t) = ^{e~^[l -(6 + j^—^ ) cothe + coth 2 e] 
-2\t u 1 o i n ^ 2 — At 2 cos(2fit) 



+e~ 2M coth e [(5 + — - 2 coth e 



5(1 -r 2 ) ' VI 2 

1 ^sin(2ta) 2r/i^(l-cos(2^)) 1 2 

" 5ii^) ] ^r- + nwi-r 2 ] + coth e} - (19) 
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In the limit of long times Eq. (fl9l) yields 



so that we obtain 



which for high T becomes 



h 2 

a(oo) = — coth 2 e, (20) 
Tioj 

5qd(oo) = tanh ^p ( 21 ) 

<W°°) = (22) 
To illustrate the dependence on temperature and time of the degree of QD, we represent 
it in Fig. 1. We see that in general Sqd decreases, and therefore QD becomes stronger, 
with time and temperature, i.e. the density matrix becomes more and more diagonal at 
higher T and the contributions of the off-diagonal elements get smaller and smaller. At 
the same time the degree of purity decreases and the degree of mixedness increases with 
T. For T = the asymptotic (final) state is pure and Sqd reaches its initial maximum 
value 1. 5qd = when the quantum coherence is completely lost, and if 5qd = 1 there 
is no QD. Only if 5qd < 1 we can say that the considered system interacting with the 
thermal bath manifests QD, when the magnitude of the elements of the density matrix 
in the position basis are peaked preferentially along the diagonal x = x'. In Fig. 2 we 
represent the density matrix in coordinate representation ffTSl) at the initial and final 
moments of time. The values of density matrix along the diagonal x = x' represent the 
probability of finding the system in this position, while off-diagonal values represent 
the correlations in the density matrix between the points x and x'. For simplicity, in 
Fig. 2 we consider zero values for the initial expectations values of coordinate and 
momentum, so that the density matrix is centered in origin. Dissipation promotes 
quantum coherences, whereas fluctuation (diffusion) reduces coherences and promotes 
QD. The balance of dissipation and fluctuation determines the final equilibrium value 
of 5 QD [2Q\. 

In order to obtain the expression of the decoherence time, we consider the coeffi- 
cient 7 (TTT1) . which measures the contribution of non-diagonal elements in the density 
matrix f)16p . For short times (At <C 1, Qt <C 1), we have: 

2 2 
TTLUJ V V LOT 

7(*) = -^{1 + 2[A(* + cothe + „{S - j^—^) cothe - X — fx — ^_]0.(23) 

From here we deduce that the quantum coherences in the density matrix decay expo- 
nentially and the decoherence time scale is given by 

2 2 

tdeco = {2[A(<5 + T71 —) cothe + ft(6 - —) cothe - A - fj, - - R - ]}~ 1 .{2A) 



5(1 -r 2 ) 7 ^ v 5(1 -r 2 ) 7 r 5^/T 



r 
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Figure 1: Degree of quantum decoherence 5qd versus temperature T (through C = 
coth (hco/2kT)) and time t for A = 0.2, y. = 0.1, S = 4, r = 0. 




Figure 2: Density matrix p in coordinate representation for A = 0.2, y — 0.1, 8 — 4, r = 
: (a) |p| at the initial time t = 0; (b) for C = 10. 
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The decoherence time depends on the temperature T and the coupling A (dissipation 
coefficient) between the system and environment, the squeezing parameter 5 and the 
initial correlation coefficient r. We notice that the decoherence time is decreasing with 
increasing dissipation, temperature and squeezing. 
For r = we obtain: 

tdec ° = 2(A + ^)(5cothe-l) (25) 
and at temperature T = [here we have to take \i = due to Eq. (JSJ)], this becomes 

1 

tdeco ~ 2A(5-1)" (26) 

We see that when the initial state is the usual coherent state (5 = 1), then the deco- 
herence time tends to infinity. This corresponds to the fact that for T = and 5=1 
the coefficient 7 is constant in time, so that the decoherence process does not occur in 
this case. 

At high temperature, expression (124p becomes 

r 2 r 2 IkT 

= {2|A(i + W^ s ) ) + *< 4 - W^) ] '"ST }_V (27) 

If, in addition r = 0, then we obtain 

tdeco = A(X + fj)6kT' (28) 

The generalized uncertainty function a(t) ffl9|) has the following behaviour for 
short times: 

a(t) = \{1 + 2[\{6 + gfi^x ) cothe + p{6 - ^^y) cothe - 2\}t}. (29) 

This expression shows explicitly the contribution for small time of uncertainty that is 
intrinsic to quantum mechanics, expressed through the Heisenberg uncertainty prin- 
ciple and uncertainty due to the coupling to the thermal environment. From Eq. 
( l29l) we can determine the time td when thermal fluctuations become comparable with 
quantum fluctuations. At high temperature we obtain 

1 1 9kT 

'-^(^sa^'+^-^a^)]^" 1 - < 30) 

As expected, the decoherence time tdeco has the same scale as the time td after which 
thermal fluctuations become comparable with quantum fluctuations [201 I2T] . The 
values of t deC o and td become closer with increasing temperature and squeezing. 
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When t 3> t re i, where t re \ ~ A -1 is the relaxation time, which governs the rate 
of energy dissipation, the particle reaches equilibrium with the environment. In the 
macroscopic domain QD occurs very much faster than relaxation, so that for all macro- 
scopic bodies the dissipation term becomes important much later after the decoher- 
ence term has already dominated and diminished the off-diagonal terms of the density 
matrix. We remark also that td eC o can be of the order of t re \ for sufficiently low tem- 
peratures and small wave packet spread (small squeezing coefficient). 

4 Equations of motion for two independent har- 
monic oscillators 

We are now interested in the generation of entanglement between two harmonic os- 
cillators due to the back-action of the environment on the subsystem. Since the two 
harmonic oscillators interact with a common environment, there will be induced cou- 
pling between the two oscillators even when initially they are uncoupled. Thus, the 
master equation for the two harmonic oscillators must account for their mutual in- 
teraction by their coupling to the environment. We shall study the dynamics of the 
subsystem composed of two identical non-interacting (independent) oscillators in weak 
interaction with a large environment, so that their reduced time evolution can be de- 
scribed by a Markovian, completely positive quantum dynamical semigroup, like in 
the previous case of one harmonic oscillator. 

If $ t is the dynamical semigroup describing the time evolution of the open quan- 
tum system in the Heisenberg picture ($ 4 is the dual of the dynamical semigroup $ t 
which describes the time evolution in the Schrodinger picture given by Eq. (CQ)), then 
the master equation is given for an operator A as follows [H] [12] : 

= l -[H, $ t (A)) + 1 J2(V^ t (A), V S \ + [Vl MAM). (31) 

Here, again H denotes the Hamiltonian of the open quantum system and Vj, Vj which 
are operators defined on the Hilbert space of H, represent the interaction of the open 
system with the environment. Being interested in the set of Gaussian states, we 
introduce those quantum dynamical semigroups that preserve that set. Therefore H is 
taken to be a polynomial of second degree in the coordinates x, y and momenta p x ,P y 
of the two quantum oscillators and Vj, Vj are taken polynomials of only first degree in 
these canonical observables. Then in the linear space spanned by the coordinates and 



10 



momenta there exist only four linearly independent operators ^-=1,2,3,4 [22J: 

Vj = a xj p x + a yj p y + b xj x + b yj y, (32) 



where a X j, a y j, b X j } b y j e C and 



Vj = a ljVx + a*jp y + b* xj x + b yj y 



(33) 



where * denotes the complex conjugation. The Hamiltonian H of the two uncoupled 
identical harmonic oscillators of mass m and frequency u is given by 



H 



(34) 



The fact that <& t is a dynamical semigroup implies the positivity of the following 
matrix formed by the scalar products of the four vectors a^, & y , b x , b y whose entries 
are the components a x j, a y j, b X j, b y j, respectively: 



h 

2 



/ ( a x a x) (&xb x ) ( a x a j/) ( a a;by) \ 

(b^a*) (b x b x ) (b x ay) (b x b v ) 



[ a y a x ) 



{a y h x ) 



\ a y a y) 



(ayby) 



(35) 



\(b y a x ) (b y b x ) (byOy) (b y b y ) ) 

For simplicity we take this matrix of the following form, where all coefficients D xx , D xpx ,... 
and A are real quantities: 



/ D xx 
-D xpx + ihX/2 



-A 



V 



D xy 



xpx -ihX/2 
Dp x p x 



ypx 



D 

-D 
D 



xy 

ypx 



Dr 



D xpy 

Dp xP y 

-D Wy - ihX/2 
D 



(36) 



'xp y ^PxPy ~Dyp y + ihX/2 ^PyPy 

It follows that the principal minors of this matrix are positive or zero. From the 
Cauchy-Schwarz inequality the following relations for the coefficients defined in Eq. 
(I36p hold (from now on we put, for simplicity, h = 1): 

A 2 

-Ti D xx D PyPy — L xp 



D xx D yy 



D% > 0, D xx D PxPx - D\^ > ^, D xx D PvPv - Dt > 0, 



'xx-^p x p x 



) 2 

XPx 



D yy D PxPx 



D z ypx > 0, D yy D PyPy 



A 2 



D z > — D D — D > 

yPy — 4 ' PxPx 'PvPy ^Px'Pv — 



(37) 
(38) 



The matrix of the coefficients (I3"6"|) can be conveniently written as 

C\ C3 

c 3 t c 2 



(39) 



in terms of 2 x 2 matrices C\ = Cy, C2 = C2 and C3. This decomposition has 
a direct physical interpretation: the elements containing the diagonal contributions 
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G\ and C*2 represent diffusion and dissipation coefficients corresponding to the first, 
respectively the second, system in absence of the other, while the elements in C3 
represent environment generated couplings between the two, initially independent, 
oscillators. 

The variance and covariance of self-adjoint operators A\ and A 2 can be written 
with the density operator p, describing the initial state of the quantum system, as 
follows: 



*AiA a (t) = -Tr(p$ t (AiA 2 + A 2 A 1 )). 
We introduce the following 4x4 covariance matrix: 



By direct calculation we obtain |22j : 

— = Yo + oY T + 2D, 
dt 



( 


@XX 


®xpx 


@xy 


a xp y 


\ 




^Xpx 


&PXPX 


a ypx 


OpxPy 




\ 


@xy 


a ypx 


a yy 


a yp v 




O'xpy 


a PxPy 


a yp y 


°PyPy 


) 



where 



/ -A 



Y 



V 



-mu 





1/m 
-X 












—muj 






1/m 
-A J 



D is the matrix of the diffusion coefficients 



D 



^ Dxx 


D X p x 


D X y 


D X p y 


\ 


Dxp x 


Dp x p x 


Dypx 


DpxPy 






Dyp x 


Dyy 


Dypy 




\D xPy 


DpxPy 


Dijpy 


DpyPy 


j 



(40) 



(41) 



(42) 



(43) 



(44) 



and Y T is the transposed matrix of Y. The time-dependent solution of Eq. 
given by [22] 

cr(t) = M(t)(a(0) - a(oo))M T (t) + a (00), 



is 



(45) 



where M(t) = exp(tY). The matrix M(t) has to fulfil the condition lim^oo M(t) = 0. 
In order that this limit exists, Y must only have eigenvalues with negative real parts. 
The values at infinity are obtained from the equation 



Fa(oo) + a(oo)y T = -2D. 



(46) 
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5 Environment induced entanglement 



The two-mode Gaussian state is entirely specified by its covariance matrix a (JJT 
which is a real, symmetric and positive matrix with the following block structure: 



a 



A C 
C T B 



(47) 



where A, B and C are 2x2 matrices. Their entries are correlations of the canonical 
operators x,y,p x ,p y , A and B denote the symmetric covariance matrices for the in- 
dividual reduced one-mode states, while the matrix C contains the cross- correlations 
between modes. The entries of the covariance matrix depend on Y and D and can 
be calculated from Eq. (1461) . To simplify further the calculations, we shall consider 
environments for which the two diagonal submatrices in Eq. fl39|) are equal: C\ — C2, 



so that D xx = Dyy,D xpx = D yPy ,D PxPx = D PyPy . In addition, in the matrix C 3 we 
take D xpy = D yPx . Then both unimodal covariance matrices are equal, A — B and 
the entanglement matrix C is symmetric. With the chosen coefficients, we obtain the 
following elements of the asymptotic entanglement matrix C: 



n ( s m 2 (2A 2 + u 2 )D xy + 2m\D xPy + D PxPy 

^ (00) = 2m 2 A(A 2 + ^ 2 ) ' (48) 



/ n / \ -m 2 uj 2 D X y + 2mXD xp + D PxP 
a xpy (oo) = a yPx (oo) = 2m(X 2 + J) " ' (49) 



, , m 2 u 4 D xy - 2mu 2 \D xpy + (2A 2 + u 2 )D PxPy 
°™ (00) = 2A(A 2 + ^ 2 ) (50) 



and of the matrices A and B: 

( v , v m 2 {2\ 2 + uj 2 )D xx + 2m\D xpx + D PxP , 
a xx (oo) = a yy (oo) = —^—^ ~ ' (51) 



/ \ , N -m 2 u 2 D xx + 2m\D xpx + D PxPx 
a xPx {oo) = a yPy {oo) = 2m(A 2 + ^ ~ > ( 52 ) 



2 4 

m 00 



D xx - 2muj 2 \D xpx + (2A 2 + oo 2 )D p 



<W«(°°) = °PyPyM = 2A(A 2 + u; 2 ) ' ( ' 

With these quantities we calculate the determinant of the entanglement matrix: 

d6t ° = A\\^+^) [{mU}2Dxy + ^ Dp ^ )2 + AX2 ( D *v D P*Py ~ D %y)\- ( 54 ) 



13 



It is very interesting that the general theory of open quantum systems allows couplings 
via the environment between uncoupled oscillators. According to the definitions of the 
environment parameters, the diffusion coefficients above can be different from zero and 
can simulate an interaction between the uncoupled oscillators. Indeed, the Gaussian 
states with det C > are separable states, but for det C < 0, it may be possible that 
the asymptotic equilibrium states are entangled, as it will be shown in the following. 

On general grounds, one expects that the effects of decoherence, counteracting 
entanglement production, be dominant in the long-time regime, so that no quantum 
correlation (entanglement) is expected to be left at infinity. Nevertheless, there are 
situations in which the environment allows the presence of entangled asymptotic equi- 
librium states. In order to investigate whether an external environment can actually 
entangle the two independent systems, we can use the partial transposition criterion 
P, El]: a state results entangled if and only if the operation of partial transposi- 
tion does not preserve its positivity. Simon [10] obtained the following necessary and 
sufficient criterion for separability: S > 0, where 

S = det A det B + (- - | det C\f - Tr[AJCJBJC r J] - -(det A + det B) (55) 
and J is the 2x2 symplectic matrix 

J=(- iJ)- < 56 > 

In order to analyze the possible persistence of the environment induced entangle- 
ment in the asymptotic long-time regime, we consider the environment characterized 
by the following values of its parameters: m 2 uj 2 D xx = D PxPx , D xpx = 0, m 2 uj 2 D xy = 
D PxP . In this case the Simon expression fl55l) takes the form: 

(m 2 u 2 (D 2 -D 2 ) D 2 l\ 2 m 2 u 2 D 2 xx D 2 / s 

\ A 2 X 2 + u 2 A J \ 2 (\ 2 + lu 2 ) 

For environments characterized by such coefficients that the expression (1571) is negative, 
the asymptotic final state is entangled. In particular, if D xy = 0, we obtain that S < 0, 
i.e. the asymptotic final state is entangled, for the following range of values of the 
coefficient D xpy characterizing the environment: 

mujD xx 1 D xpy muD xx 1 

2 < 7FT^ < ^\ _ + 2' (58) 

where the coefficient D xx satisfies the condition muoD xx /\ > 1/2, equivalent with the 
unimodal uncertainty relation. If the coefficients do not fulfil the inequalities (1581) . 
then S > and therefore the asymptotic final state of the considered bipartite system 
is separable. 
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6 Summary 



We have studied QD with the Markovian equation of Lindblad for a system consist- 
ing of an one-dimensional harmonic oscillator in interaction with a thermal bath in 
the framework of the theory of open quantum systems based on completely positive 
quantum dynamical semigroups. In the same framework we investigated the existence 
of the asymptotic quantum entanglement for a subsystem composed of two uncoupled 
identical harmonic oscillators interacting with an environment. 

(1) We have shown that QD in general increases with time and temperature. For 
large temperatures, QD is strong and the degree of mixedness is high, while for zero 
temperature the asymptotic final state is pure. With increasing squeezing parameter 
and initial correlation, QD becomes stronger, but the asymptotic value of the degree of 
QD does not depend on the initial squeezing and correlation, it depends on temperature 
only. QD is expressed by the loss of quantum coherences in the case of a thermal bath 
at finite temperature. 

(2) We determined the general expression of the decoherence time, which shows 
that it is decreasing with increasing dissipation, temperature and squeezing. We have 
also shown that the decoherence time has the same scale as the time after which 
thermal fluctuations become comparable with quantum fluctuations and the values of 
these scales become closer with increasing temperature and squeezing. 

(3) By using the Peres-Simon necessary and sufficient condition for separability 
of two-mode Gaussian states, we have shown that for certain classes of environments 
the initial state evolves asymptotically to an equilibrium state which is entangled, 
i.e. there exist non-local quantum correlations for the bipartite states of the two 
harmonic oscillator subsystem, while for other values of the coefficients describing the 
environment, the asymptotic state is separable. 

The obtained results can represent a useful basis for the description of the connec- 
tion between uncertainty, decoherence and correlations (entanglement) of open quan- 
tum systems with their environment. Due to the increasing interest manifested to- 
wards the continuous variables approach [23] to the theory of quantum information, 
these results, in particular the possibility of maintaining a bipartite entanglement in 
a diffusive-dissipative environment even for asymptotic long times, could be useful for 
both phenomenological and experimental applications in the field of quantum infor- 
mation processing and communication. 
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